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l/^ . Abstract 
CN ■ 

The publicly available Super-Kamiokande 10 day plot data have 
■ been searched for a periodic modulation in the range between 10 and 

500 days, with the periodogram method. A peak is observed for T = 
13.75 days, significant at 98.9 % C.L. I discuss possible ways to check 



^ • this result and, if confirmed, to investigate its origin. 
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^ : 1 Introduction 

In recent years, we have witnessed a dramatic improvement in our assessment 
of the solar neutrino problem. The very high statistics of solar neutrino interac- 
tions in ordinary water, obtained in the Super-Kamiokande experiment lUl^l, 
together with results from the chlorine gallium j31 El E] and heavy wa- 
ter 13 IH] detectors, have demonstrated that neutrinos undergo flavor conversion 
between the solar core and the earthly experiments. In parallel and at about 
the same time, the enormous progress in helioseismology and in solar models 
led to a firmer and refined estimate of the emitted solar neutrino flux [U] , that 
agrees perfectly with the measurements. A few domains were singled out in 



the Am? vs. tanO parameter plane (see e.g. pOj), with the so-called Large 
Mixing Angle solution being the favored one. The very recent observation of 
reactor neutrino oscillations by the KamLAND group jTT] has further pinned 
down this solution, although the allowed domain in parameter space is now 
comprised of two disjoint sub-domains (see e.g. |12|). 

In the present standard view, no large periodic modulation is expected in 
the solar neutrino signal. Two small effects might be observed, with periods 
of one year and one day. The obvious yearly modulation is due to the Earth's 
orbit excentricity, e = 0.0167, which induces a flux variation 4e/(l — e)^ ~ 7% 
between aphelion and perihelion. There is now some evidence for this flux 
modulation, at a modest significance level |T|. The possible daily periodicity 
stems from neutrino oscillations inside the Earth :13^. The amplitude of this 
so-called day/night effect depends strongly on the neutrino masses and mix- 
ings. The day/night comparisons by the Super-Kamiokande 0121 and SNO jH] 
collaborations helped to reject some solutions to the solar neutrino problem, 
and to narrow down the allowed parameter space domain. 

Searches for other periodicities have also been performed, such as correla- 
tions with the 11-year solar cycle {e.g. Jl]), or the solar synodic period (see 
e.g. |15| ITB] and references therein); it seems fair to say that no consensus 
exists on the presence of such periodic modulations. 

Very recently, in December 2002, the Super-Kamiokande collaboration has 
released its solar neutrino flux data in 10 and 45 day bins over the period May 
1996 to July 2001, corresponding to 22,400 solar neutrino interactions (see 
http://www-sk.icrr.u-tokyo.ac.jp/sk/lowe/index.html). This makes it 
possible to perform a systematic search for any periodicity in this data set, 
not restricted to the four periods listed above. The present paper reports the 
results of such a search in the period range from 10 to 500 days, as well as 
evidence for a periodic modulation, at T = 13.755 ± 0.011 days. 1 am not 
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aware of any previous systematic periodicity search on solar neutrino fluxes in 
a wide period range, and especially below 18 days^ 

2 The data sample and the search for period- 
icity 

I have used the Super-Kamiokande solar neutrino flux in 10 day bins only, 
in order to probe the shorter periods. There are = 184 such data points, 
for each of which are provided the date and time of the beginning and end of 
the 10-day bin (in Japan Standard Time JST), the value of (f, the average 
squared Sun-Earth distance in this bin, and the boron-8 solar neutrino flux, 
together with errors, uncorrected for this distance variation. The relative flux 
errors are typically between 10 and 20%. 

The total time span of the data is 1,864 days (184 bins of 10.1 days in 
average), corresponding to an actual live time of 1,496 days. The bin size is 
rather constant, as only 13 bins are shorter than 8 or longer than 12 days, 
and the standard deviation of the bin size distribution is 1.27 day. The actual 
distribution of dead time is not provided, so that I have chosen to assign the 
neutrino flux value to the center of each bin. I have then corrected this flux 
value for the Sun-Earth distance, using in effect a flux rescaled at 1 AU. (I 
have used the d"^ factor provided in the data table, and have checked that it 
agrees with that computed using the JPL DE406 Solar System ephemerides 
available at |http : / / ssd • jp l ■ nasa . gov/ . ) Finally, I subtract in each bin the 
average flux value given in the same table. 

I have then performed a search for periodicity on these "flux excursions" , 

^Half a day after this preprint was first posted, I became aware of a similar search 
by P.A. Sturrock and D.O. CaldweU, through an AAS meeting abstract, posted at 
[http: //www. aas . org/publications/baas/v34n4/aas20 1/64 7. htm . I very much thank 
T. Lasserre for pointing this out to me. The periodicity found in the present paper agrees 
perfectly with that given in this conference abstract. 
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using a slight variant of the periodogram method f7|^]. All flux values were 
used with equal weights, i. e. ignoring their errors. The frequency range from 
0.002 to 0.100 d^^ was searched, with a frequency step of 0.0001 d^^; this 
corresponds to a period range between 500 and 10 days. (The most significant 
peaks were then scanned again, with a 10 times smaller step, 0.00001 d~^.) 
The resulting power spectrum is shown in Fig. 1, where the plotted variable is 

where is the flux measured at time ti, the average flux and ly is the 
probed frequency; the sums run over the data points. 

Only one prominent peak is found, at a frequency u = 0.07270 d^^ ~ 
842 nHz, corresponding to T = 13.75 ± 0.05 days. (This error is only an 
estimate, from the width of the periodogram peak.) The height of the peak 
is Zq = 10.4, corresponding to a probability p = e~^° ~ 3 ■ 10^^. The z 
distribution is shown in Fig. 2. In order to estimate the probability that 
such a peak be compatible with a statistical fluctuation, one has to know the 
number of independent frequencies tested by the periodogram analysis. For 
the frequency range tested here, essentially between and twice the Nyquist 
frequency, this number is close to 2A^, i.e. twice the number of data points 
(see e.g. ^H]; section 13.8, and references therein). The probability to obtain a 
peak height greater than zq from 2A^ tested frequencies is 1 — (1— p)^^ ~ 2Np ~ 
1.1%. Thus the observed peak corresponds, at 98.9% C.L., to a periodicity in 
the Super-Kamiokande 10 day plot data (as of December 2002). The second 
highest peak of the periodogram, with zi = 6.4 , is perfectly compatible with 
statistical noise at 47% CL. 

Note that the sampling, though quite regular, is sufficiently variable that 
no aliasing is observed in the periodogram. 

The folded flux values are plotted versus their phase in Fig. 3. The phase 
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here is defined as the fraction of a period T counted from an arbitrary origin 
(Modified Juhan Day 0, or JD 2,450,000 = 1995, Oct 9.50 UTC); it is a number 
between and 1. The modulation is clearly visible on Fig. 3. 

In order to better estimate the parameters and their errors, I have fitted 
the 184 data points cr^.) with a sinusoidal fiux modulation, 

^ r. . • ti- 1200 , 

$0 [1 + Asin{2Ti^^ + 0i2oo)], 

where $o is the time-averaged neutrino fiux, A the relative signal mod- 
ulation, T the period and 0i2oo the phase at a reference date, MJD 1200.0, 
chosen close to the median time of the SK data. In this fit, the points were 
weighted using the inverse squared errors quoted by the Super-Kamiokande 
collaboration. The resulting is 165.8 for (184 — 4) degrees of freedom. 
The fitted parameters values are : $o = 2.329 ± 0.02410^ ^B-v 01x1-"^ s''^ , 
A = 0.067 ± 0.014, T = 13.755 ± 0.011 days and 0i2oo = 0.55 ± 0.22 rad. (In 
comparison, the of a constant fit is 187.7 for 183 d.o.f.) 

It may seem strange to be able to detect a period that is only 35 % longer 
than the average time bin size. The neutrino flux in each bin is integrated over 
the bin width, so that any periodic signal will be washed out when integrated 
over ~ 3/4 of a period. In the case of a purely sinusoidal modulation, with 
amplitude A relative to the average fiux, the measured amplitude after signal 
averaging over a random 3/4 fraction of a period is A! — (2^/^/37r)A ~ 0.30>1. 
I will discuss this further in the next section. 

Note that, for such an integrated signal, the periodogram search has a 
much reduced sensitivity to periods close to the time sampling, 10.1 days, 
and twice the time sampling, 20.2 days. Thus the number of independent 
frequencies, 2A'", that was used above is probably a bit overestimated, and the 
signal significance under-estimated. A full Monte-Carlo simulation would be 
necessary to quantify this. 
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In order to search for a possible bias, either on the 13.75 day period or 
on its significance level, I have performed a Monte-Carlo simulation gener- 
ating 10,000 experiments with exactly the same 10 day time bins as Super- 
Kamiokande, and flux values taken at random from a normal law with the 
same average (2.35 ■ 10^ z/cm~^ s~^) and the same dispersion (0.33 ■ 10^) as the 
184 time bins. The observed z distribution is in perfect agreement with an ex- 
ponential law, and no significant bias is observed near 13.75 days, or anywhere 
else in the periodogram. Out of the 10,000 simulated data sets, 123 have a 
maximum z value larger than = 10.4 . This agrees with our confidence level 
estimate of 98.9%. 

I have also performed a periodogram search on the combined results of 
the GALLEX ;4J and GNO [F collaborations over the period 1991-2001. No 
significant peak was found. This gives no information on a possible 13.75 day 
periodicity however : out of the 100 GALLEX and GNO runs, 65 correspond to 
an exposure time very close to 27.5 days, i.e twice 13.75 days, so that a much 
reduced sensitivity is expected. I have not analysed the chlorine data, as the 
time between extractions is even longer for this experiment. 

3 Discussion and Future Work 

Provided the signal presented in the previous section is not simply due to a 
statistical fluctuation, there exist many tests that could be performed to verify 
its existence. An identical analysis using the same Super-Kamiokande data, 
but binned in 3-4 days bins would already prove very useful. The statistics in 
each bin would decrease by a factor 3, and the statistical errors would increase 
typically by -\/3- The distribution of measured fluxes would hence be wider. 
But a 13.75 day periodic signal would then, in these bins, be integrated over 
1/4 of a period instead of 3/4. The amplitude of the signal would increase 
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by a factor 3, to about 20%, thus more than compensating the increase in 
the width of the distribution; the significance of any real modulation could be 
clearly established. 

It would then be interesting to check what the actual period is, as 13.75 days 
could very well be the (dominant) second harmonic of a 27.5 day period. (The 
periodogram method is not well suited in such the AoV method [20] 

could be used instead.) One should also check for aliases with periods shorter 
than 10 days. 

A search for periodicity is always more efficient starting from unbinned 
data. The ultimate verification would be to search for periodicity using infor- 
mation on each "candidate solar neutrino", i.e. the time of detection, the angle 
with respect to the Sun {cos9q) and the energy. For the Super-Kamiokande 
data, one can only speak of "candidate solar neutrinos" as the boron-8 solar 
neutrino fiux is measured as an excess over the roughly constant background 
in the cos9q distribution (see for example figure 1 of [211 )• One could use 
a modified version of a periodogram analysis, where each "candidate" would 
receive a cos6'0-dependent weight, equal to the fraction of solar neutrinos at 
its measured cosOq. (The closer 6'© to zero, the larger the weight.) 

A simplified Monte-Carlo simulation of such a weighted periodogram anal- 
ysis has been performed [22] • It shows that a search for an unknown periodic 
modulation of the fiux in the period range from a few hours to a year would 
allow to detect a 10 % modulation at better than 99 % CL. (The confidence 
level increases very fast with the amplitude.) This means that, for a fixed 
confidence level, the searched frequency interval would then be about a hun- 
dred times greater than in the present paper. Equivalently, examining in such 
an analysis only periods greater than 10 days would lead to confidence levels 
100 times better. 

In case the signal is confirmed, this would not tell us whether it is due to 
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an instrumental effect, to an actual modulation of the boron-8 neutrino flux, 
or if it is induced by a modulation in the background. The latter hypothesis 
could be tested, for example, by a periodogram analysis of events that en- 
ter the cosOq plot at negative values (the anti-solar hemisphere) : they are 
almost all background events. It would also be interesting to search for the 
same periodicity in other Super-Kamiokande data sets, such as atmospheric 
neutrinos. 

One cannot help noticing that 13.75 days is very close to half the solar 
synodic rotation period. This may turn out to be a coincidence, or not. There 
exist many solar observables that display a 13-14 day periodicity, although the 
peak in their Fourier spectrum is usually quite wide. Examples of these are 
properties of the chromosphere, solar wind characteristics at the Earth, and 
even induced geomagnetic activity (see e.g. and references therein). (The 
latter two might influence the background level in underground experiments; a 
quick bibliographic search did not allow me to find a reference in the litterature 
to a periodicity in the flux of cosmic rays.) Note however that the observed 
13.5 day periodicities have not been seen to persist over a very long time 
span, such as the 5 years of data analysed here. Moreover, the intensity of 
these modulations vary during the 11-year solar cycle, and often show phase 
changes between consecutive periods of 13.5 day variability. For that reason, I 
have tried to split the Super-Kamiokande data into various time bins to search 
for such phase shifts, with no positive results. It is likely however that the 10- 
day binned data sensitivity is not large enough to detect such effects, and that 
the periodogram is not the proper tool here. 

In case there is a connection between the observed signal and solar prop- 
erties, it would certainly be interesting to search for a correlation with the 
position of the Earth relative to the Sun's heliographic or heliomagnetic lati- 
tude. 
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Regarding other experiments, we have aheady mentioned that the galhum 
and chlorine data have a lower sensitivity, due to the long time between ex- 
tractions, and to an unfortunate coincidence. If the Super-Kamiokande signal 
is confirmed by further analyses, the matter should of course be reconsidered, 
as one would then just have to study the amplitude of a modulation at a 
known period. In contrast, all other solar neutrino experiments with "real 
time detection" such as SNO, or the forthcoming Borexino [23^, can use the 
techniques described above or similar ones to search for a signal. One would 
expect the SNO data to display the same kind of variation as the SK data, at 
the level of 20%. 

4 Summary 

I have presented the results of a search for a periodic modulation in the Super- 
Kamiokande 10 day plot data (as of December 2002). At 98.9 % C.L., a 
T = 13.75 ± 0.01 day modulation is found. This is close to half the Sun's 
synodic rotation period. I have discussed methods to soon check this result, 
and have sketched possible future directions, in case it is confirmed. 
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Figure Captions 

1. The power spectrum in the Super-Kamiokandc 10 day plot data, from 
a periodogram analysis. The levels corresponding to a 1, 5 and 20 percent 
probability of a statistical fluctuation correspond here to 2; = 10.5, 8.9 and 
7.4, respectively. 

2. The z distribution corresponding to Fig. 1. For data with no significant 
periodicity, this distribution should be exponential (oc e^^). The three z values 
in excess of 6.5 all belong to the same periodogram peak v — 0.07270 d~^. 

3. The flux values versus their phase in a period T — 13.75 d. The phase 
here is the fraction of the period T, counted from an arbitrary origin (MJD 
0.0 = JD 2,450,000). It is a number between and 1, but each data bin has 
been plotted twice, at its phase and phase + 1. 
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